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ABSTRACT 

We have used spectra obtained with X-shooter, the triple arm optical-infrared spectrograph re- 
cently commissioned on the Very Large Telescope (VLT) of the European Southern Observa- 
tory (ESO), to confirm the gravitational lens nature of the CASSOWARY candidate CSWA 20. 
This system consists of a luminous red galaxy at redshift Zabs = 0.741, with a very high veloc- 
ity dispersion, cricns — 500km s^^, which lenses a blue star-forming galaxy at z^m = 1.433 
into four images with mean separation of ~ 6". The source shares many of its properties with 
those of UV-selected galaxies at z — 2-3: it is forming stars at a rate SFR ~ 25Mq yr^^, has 
a metallicity of 1/4 solar, and shows nebular emission from two components separated by 
0.4" (in the image plane), possibly indicating a merger. It appears that foreground interstellar 
material within the galaxy has been evacuated from the sight-line along which we observe the 
starburst, giving an unextinguished view of its stars and H II regions. CSWA 20, with its mas- 
sive lensing galaxy producing a high magnification of an intrinsically luminous background 
galaxy, is a promising target for future studies at a variety of wavelengths. 

Key words: gravitational lensing - galaxies: evolution - galaxies: structure. 



1 INTRODUCTION 

The large cosmic volume surveyed by the Sloan Digital Sky Sur- 
vey (SDSS) has sparked, among many other projects, several sys- 
tematic searches for strong gravitational lens systems (Bolton et al. 
2006; Willis et al. 2006; Estrada et al. 2007; Ofek et al. 2008; Shin 
et al. 2008; Belokurov et al. 2009; Kubo et al. 2009; Lin et al. 2009; 
Wen et al. 2009). All of these studies share the dual motivation of, 
on the one hand, probing the high-mass end of the galaxy mass 
function and the underlying distribution of dark matter in the lens- 
ing galaxies and, on the other, identifying highly magnified high 
redshift sources. The latter can bring within reach of current as- 
tronomical instrumentation detailed studies of stellar populations 
and interstellar gas at high redshifts which would otherwise have 
to wait until the advent of the next generation of 30+ m optical- 
infrared telescopes (e.g. Pettini et al. 2000, 2002; Teplitz et al. 
2000; Lemoine-Busserolle et al. 2003; Small et al. 2007; Cabanac, 
Valls-Gabaud, & Lidman 2008; Siana et al. 2008, 2009; Finkelstein 



* Based on public data from the X-shooter commissioning observations 
collected at the European Southern Observatory VLT/Melipal telescope, 
Paranal, Chile. 



et al. 2009; Mainline et al. 2009; Quider et al. 2009, 2010; Yuan & 
Kewley 2009 and references therein). 

The CAmbridge Sloan Survey Of Wide ARcs in the skY 
(CASSOWARY) targets multiple, blue companions around mas- 
sive ellipticals in the SDSS photometric catalogue as likely can- 
didates for wide-separation gravitational lens systems. A com- 
prehensive description of the search strategy is given by Be- 
lokurov et al. (2009). Of the twenty highest priority CASSOWARY 
candidates, eight have so far been confirmed as gravitational 
lenses and the redshifts of both lens and source measured — see 
|http;//www.ast.cam.ac.iik/research/cassowary71 for further details. 

In this paper, we report observations of a ninth system, CAS- 
SOWARY 20 or CSWA 20 for short. As can be seen from Figure[I] 
CSWA 20 consists of four blue images around a red galaxy, in a 
configuration reminiscent of the Einstein Cross (Adam et al. 1989), 
but with a factor of ~ 3 larger separation between the images. 
The observations reported here, obtained during the commission- 
ing of the triple arm spectrograph X-shooter on the VLT, confirm 
the gravitational lens nature of the system by showing that three of 
the four blue images are of the same source, a star-forming galaxy 
at redshift Zcm = 1.433, and that the lens is a massive elliptical 
galaxy at z^bs = 0.741 (the fourth blue image was not observed. 
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Table 1. SDSS Positions and Magnitudes of Components of CSWA20 Lens System 

Image RA (J2000) Dec (J2000) u g r i z 

Lens 14 41 49.16 +14 41 20.6 24.6 ±1.1 25.1 ±0.5 22.7 ± 0.2 20.70 ±0.06 20.2 ± 0.2 

11 1441 49.38 +1441 21.3 21.46±0.15 21.60±0.05 21.54±0.07 21.87±0.17 21.76±0.63 

12 14 41 49.24 +14 41 22.8 21.89 ±0.26 21.65 ±0.06 21.40 ± 0.08 21.27 ±0.13 21.89 ±0.88 

13 14 41 48.93 +14 41 22.2 21.56 ±0.23 21.52 ±0.07 21.46 ±0.11 21.27 ±0.16 21.18 ±0.67 

14 1441 49.18 +1441 18.7 21.58±0.23 21.57±0.07 21.16±0.08 21.19±0.16 20.59 ± 0.38 



Note: Positional errors are ^ 0.1" for images il-i4, and ~ 0.1" for the lensing galaxy. 




Figure 1. Colour-composite (g, r, i) SDSS image of the CSWA20 lens sys- 
tem. Positions and magnitudes of the different components are listed in Ta- 
ble[T] 



in collaboration with ESO which was responsible for the final in- 
tegration and installation on the VLT. A full description of the in- 
strument is provided by D'Odorico et al. (2006). X-shooter con- 
sists of three echelle spectrographs with prism cross-dispersion, 
mounted on a common structure at the Cassegrain focus of the 
Unit Telescope 2 (Kueyen). The light beam from the telescope is 
split in the instrument by two dichroics which direct the light in 
the spectral ranges 300-550 nm and 550-1015 nm to the slits of the 
UV-B and VIS-R spectrographs respectively. The undeviated beam 
feeds the NIR spectrograph with wavelengths in the range 1025- 
2400 nm. The UV-B and VIS-R spectrographs operate at ambient 
temperature and pressure and deliver two-dimensional spectra on 
the 2048 x 4102 15 /im pixel E2V CCD and 2048 x 4096, 15 ^tm 
pixel MIT/LL CCD respectively. The NIR spectrograph is enclosed 
in a vacuum vessel and kept at a temperature of approximately 
80K by a continuous flow of liquid nitrogen. The NIR detector 
is a Teledyne substrate-removed HgCdTe Hawaii-2RG array, with 
2048 X 2048 18 /im pixels. The spectral format in the three spec- 
trographs is fixed; the final spectral resolution is determined by the 
choice of slit width with each spectrograph having its own slit se- 
lection device. 



but is most likely to be a fourth gravitationally lensed image of 
the z — 1.433 galaxy). Table [T] gives positions and magnitudes 
for the different components of the system. The four images of the 
source are all of similar magnitudes (e.g. g — 21.52-21.65) and 
all have very blue colours, with —0.14 ^ {u ~ g) ^ +0.24 and 
+0.06 ^ ((? — r) ^ +0.41. In contrast, the lens has red colours 
with [g -r) = +2.4 and (r - i) = +2.0. 

In Section 2, we give a concise description of the X-shooter 
instrument and details of the observations and data reduction. 
Section 3 describes a simple lensing model we have applied to 
CSWA20, while Sections 4 and 5 present our measurements of the 
spectra of the lens and source respectively. We summarise our re- 
sults in Section 6. Throughout this paper we use a '737' cosmology 
with Ho = 70km s"^ Mpc"\ = 0.3 and Oa = 0.7. 



2 OBSERVATIONS AND DATA REDUCTION 
2.1 X-shooter 

X-shooter is the first of the second generation VLT instruments to 
be made available to the ESO community; it was built by a consor- 
tium of institutes in Denmark, France, Italy and the Netherlands, 



2.2 Observations 

After two commissioning runs with the UV-B and VIS-R spectro- 
graphs in November 2008 and January 2009, the instrument was 
operated in its full three-arm configuration in two further commis- 
sioning runs in March 2009 and May 2009 during which the ob- 
servations of CSWA 20 were obtained. This lens candidate was se- 
lected as a good test of the instrument performance for faint galaxy 
studies. Different observation strategies were attempted, as detailed 
in Table (2] In March 2009 the 11" long entrance slit was rotated 
to a position angle on the sky PA= 117° and positioned so as 
to record simultaneously the spectra of images il and 12 (see Fig- 
ureO. In May 2009 the slit was placed across image 12 and the 
lensing galaxy (R4= 11°); this set-up also captured some of the 
light from image 14. For all observations the slit widths were 1" 
(UV-B), 0.9" (VIS-R), and 0.9" (NIR); the corresponding resolv- 
ing powers are ii = A/ A A = 5100, 8800, and 5600, sampled 
with 3.2, 3.0, and 4.0 wavelength bins respectively (after on-chip 
binning by a factor of two of the UV-B and VIS-R detectors). The 
observations used a nodding along the slit approach, with an offset 
of 4" between individual exposures. All calibration frames needed 
for data processing were obtained on the day after the observations. 
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Table 2. Details of X-shooter Observations of CSWA 20 



Date (UT) 


Exp. Time (s) 


Slit PA (°) 


Comments 


2009 March 19 


2 X 1500 (UV-B, VIS-R) 


117 


Two exposures on targets, nodded. 




2 X 1600 (NIR) 




Slit across il and 12 (see Figure[T). 


2009 May 05 


4 X 1500 (UV-B, VIS-R, NIR) 


11 


Two exposures on targets, two on sky. 








Slit across 12, lensing galaxy, and marginally 14. 



2.3 Data Reduction 

The spectra were processed with a preliminary version of the X- 
shooter data reduction pipeline (Goldoni et al. 2006). Pixels in the 
two-dimensional (2D) frames are first mapped to wavelength space 
using calibration frames. Sky emission lines are subtracted before 
any resampling using the method developed by Kelson (2003). 
The different orders are then extracted, rectified, wavelength cal- 
ibrated and merged, with a weighted average used in the overlap- 
ping regions. The final product is a one dimensional, background- 
subtracted spectrum and the corresponding error file. Intermediate 
products including the sky spectrum and individual echelle orders 
are also available. We followed all of these steps, although we used 
standard IRAF tools for extracting the ID spectra (using a prede- 
fined aperture) while this aspect of the pipeline data processing was 
being refined. 

The VIS-R and NIR spectra were corrected for telluric absorp- 
tion by dividing their spectra by that of the 08.5 star Hip 69892, 
observed with the same instrumental set-up and at approximately 
the same airmass as CSWA 20. Absolute flux calibration used as 
reference the Hubble Space Telescope white dwarf standard GD 71 
(Bohlin et al. 2001) whose spectrum was recorded during the same 
nights as CSWA 20. 



3 LENSING MODELS 

As discussed in detail in Sections |4] and |5] the X-shooter spectra 
show the lens to be an absorption line galaxy at Zabs = 0.741, 
and three of the four blue images (il, 12, and 14) to be those of an 
emission line galaxy at Zem = 1.433. With the assumption that 
13 is also a gravitationally lensed image of the same galaxy, we 
can use these redhifts together with the data in Table [T] to develop 
some simple lensing models for CSWA 20. The aims are to obtain 
estimates of the enclosed mass within the images and hence the 
velocity dispersion of the lensing galaxy, as well as estimates of 
the total magnification of the images. 

Before formal modelling, let us begin with a very simple 
model to estimate rough, order of magnitude effects. Suppose the 
lens is an isolated singular isothermal sphere with a constant veloc- 
ity dispersion a. The lensing properties of this model are discussed 
by Schneider, Ehlers & Falco (1992), who show that the typical 
deflection is: 

Ae = l'.'15f (1) 

V200kms-V ) 

where _Dds is the angular diameter distance between deflector and 
source, whilst is the distance between observer and source. A 
simple estimate of the velocity dispersion can be immediately ob- 
tained by requiring that the isothermal sphere deflection given by 
eq.lUll reproduce the observed deflections (AS ~ 3") for the dis- 
tances Dds and implied by the values of Zabs and Zcm in our 




Figure 2. A PixeLens model of CSWA 20 with image locations (red dots) 
and source locations (black dots). The left panels shows the Fermat surface, 
from which the parity of the images can be inferred. The right panel shows 
the pixellated mass distribution in units of the critical surface mass density. 



cosmology. This already suggests that the lensing galaxy is very 
massive with a ~ 500 km s^^. 

Of course, mass distributions in nature are not spherically 
symmetric. Ellipticity occurs both because of the intrinsic flatten- 
ing of the lens galaxy and because of the external tidal shear gen- 
erated by neighbouring galaxies. We therefore wish to supplement 
our simple model with more realistic, flattened mass distributions. 
These can reproduce not just the separations, but also the detailed 
positions of all four putative images. 

A powerful way of modelling gravitational lenses is to pixel- 
late the projected mass distribution of the lensing galaxies into 
tiles. Mass can be apportioned to each tile. The mass on the tiles 
is unknown, but fixed by requiring that the mass distribution repro- 
duce the images with the observed parities and locations. Of course 
the problem is then under-determined, as there are many more un- 
knowns than constraints, but can be regularised by requiring that 
the mass distribution is isothermal-like. This simple idea has been 
developed by Williams & Saha (2004) into the PixeLens code, 
which provides flattened generalisations of the isothermal sphere. 

In detail, the solution space for the masses on the tiles is sam- 
pled using a Markov chain Monte Carlo method. We typically gen- 
erate an ensemble of 1000 models that reproduce the input data, 
which in our case are the image locations given in Table [T] As the 
constraint equations are linear, averaging the ensemble also pro- 
duces a solution, which is displayed in Figure[2l 

The left panel shows the arrival time surface. The images lie 
at the stationary points of the surface and are marked with red dots. 
In the model, we see that the images il and 13 correspond to local 
minima and so are positive parity, whereas 12 and 14 correspond to 
saddles and so are negative parity. The model predicts an additional, 
de-magnified central image which - as is customary - is too faint 
for detection. The right panel shows the contours of projected mass 
of the absorption line galaxy (in terms of the critical surface mass 
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Figure 3. Three spectra are shown in this Figure. Black: Portion of the X-shooter spectrum of the lensing galaxy in CSWA 20, boxcar smoothed with a three- 
pixel wide filter. Blue: The spectrum of SDSSJ010354. 1+144814.1 the galaxy with the largest velocity dispersion (tj = 530 ± 60km s"'^) in the sample of 
Bemardi et al. (2006). Red: Model spectrum of a 12 Gyr old single burst of star formation, from Maraston et al. (2009). The inset shows a detailed view of the 
spectral region encompassing the Call H & K lines, together with our Gaussian fit. 



density). This is mildly elliptical in the inner parts, and falls off in 
an isothermal-like manner. 

From the PixeLens model, we compute the (cylindrical) 
mass within the Einstein radius (21.1 kpc) as M — (4.4 ± 0.4) x 
10^^ Mq. This implies a velocity dispersion for the lensing galaxy 
of 530 ± 20 km s ^ ^ , if it is isothermal. As we shall see in Section|4l 
such a high velocity dispersion is consistent with the widths of the 
absorption lines in the X-shooter spectrum of the lens. The errors 
have been calculated using PixeLens to generate 1000 models 
that reproduce the input data. From these distributions, we calcu- 
late the error on the enclosed mass by constructing 68 per cent con- 
fidence limits. The error on the enclosed mass gives an error on the 
velocity dispersion via the Jeans equations. The total magnification 
is rather poorly constrained with the present data, but most of the 
models have a total magnification (for the four images) of factors 
between ~ 3 and ~ 6 



4 THE LENS: A MASSIVE GALAXY AT Zabs = 0.741 

The spectrum of the deflector galaxy, reproduced in Figure[3] shows 
the signatures of an old stellar population at a redshift 2:abs = 
0.741. The galaxy has a prominent 4000 A break, and strong 
Can H & K, G-band and Mg6 AA5167.3, 5172.7, 5183.6 absorp- 
tion features. H/3 absorption is also present. The signal-to-noise ra- 
tio of the data is only modest, but Gaussian fits to the Ca II H&K 
absorption (see inset in Figure |3j result in an estimate of the ve- 
locity dispersion a = 495 ± 54 km s"'^ after correcting for instru- 
mental broadening (which is minimal, since ainstr — 15 km s^^). 
Additional observations with higher signal-to-noise ratio are re- 
quired to provide a significantly improved measure of the velocity 
dispersion. The overall spectral energy distribution at optical wave- 
lengths is similar to those of luminous red galaxies (LRGs) at lower 
redshifts, as can be appreciated from the comparison in Figure [3] 
with the SDSS spectrum of J010354.1-I-144814.1, the most massive 



galaxy found by Bemardi et al. (2006) with no obvious indication 
that its large velocity dispersion (a — 530 ± 60 km s~^) may be 
due to the superposition of two objects along the line of sight. Also 
shown in Figure[3]is the model spectrum computed by Maraston et 
al. (2009) for a 12 Gyr old single burst of star formation with solar 
metallicity; this comparison further illustrates the very red contin- 
uum (longward of the G-band) of the deflector galaxy in CSWA 20. 



Given the relatively high redshift, and associated (1 + z)* 
surface brightness dimming, it is not surprising that the galaxy 
is barely detected in the SDSS r-band image and that the i-band 
magnitude is poorly determined. Only the core of the galaxy is 
evident in the SDSS images, with an apparent effective radius of 
< 3 kpc, and deeper, higher-resolution imaging is needed to de- 
termine its structural properties. However, we can put some con- 
straints on its total luminosity from the X-shooter spectral data. The 
integrals of the deflector's spectral energy distribution multiplied 
by the SDSS r-band and i-band response curves give r — 21.7 and 
i — 20.1 respectively. Given the slit width of 0.9", a mean seeing 
of FWHM — 0.65" , and assuming an effective radius of the galaxy 
rc — 10 kpc (typical of LRGs, Bernardi et al. 2008) which corre- 
sponds to T-c = 1.1" at 2 = 0.741, we calculate that ~ 40% of 
the total galaxy light falls within the X-shooter slit. The estimate 
of the i-band magnitude is then i — 19.1. Adopting a fe-correction 
of —1.0 mag, an evolutionary correction of 0.8 mag and rest-frame 
colour (r — i) — 0.4 (appropriate to a 12 Gyr old passively evolving 
Maraston et al. 2009 model), then gives a corrected r = 18.9, and 
a corresponding absolute magnitude of Mr = —24.4. With these 
parameters, this galaxy is among the most massive and luminous 
red galaxies known (Bemardi et al. 2006). 
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Figure 4. A selection of nebular emission lines in the CSWA 20 lensed galaxy. In each panel, the black histogram shows the data, while the blue continuous 
line is the fit with the parameters given in Table|4] 



THE SOURCE: A LUMINOUS STAR-FORMING 
GALAXY AT Zem = 1.433 



As mentioned in Section 12.21 our observations of the source in 
CSWA 20 cover images il and i2 (the latter observed at two 
epochs — see Table[2l(. In order to improve the signal-to-noise ratio 
(S/N), we added together the spectra of il and i2 (after convert- 
ing the wavelengths to a vacuum heliocentric frame of reference 
and biiming to a common wavelength grid) and used this compos- 
ite spectrum in the analysis described below. The spectrum shows 
a number of narrow emission lines superposed on a weak blue con- 
tinuum (see FigureQ; the line identifications in Table[3]indicate a 
redshift Zem = 1.433. 

We also recorded image i4 on the second observing run on 
2009 May 05, but the slit position, chosen to cover image i2 and 
the Lens, only captured a small fraction of the light of i4. While 
this spectrum does show the strongest emission lines at a similar 
redshift as il and i2 (with small differences attributable to the fact 
that the image of i4 was offset relative to the slit centre, leading to 
an offset in the wavelength calibration), it was not included in the 
composite because its S/N is much lower than that of il and i2. 



5.1 Nebular Emission Lines 

From the lensing model described in Section [3] it was concluded 
that the overall magnification factor for the four images of the Ein- 
stein cross is between ~ 3 and ~ 6. It can also be seen from 
Table [T] that in the u and g filters images il and i2 account for ap- 
proximately half of the total flux, that is ii(il-l-i2) ~ M(i3-l-i4), and 



similarly for the g magnitudes. Thus, in all the following analy- 
sis we shall assume, for simplicity, that a magnification factor of 
~ 5/2 = 2.5 applies to the quantities measured for il-l-i2. 

With this assumption and the knowledge of the redshift, we 
can immediately estimate the luminosity of the source, and com- 
pare it with that of other galaxies at similar redshifts. Luminos- 
ity functions at z ^ 2 have been measured mostly in the rest- 
frame far-UV continuum, at 1700 A (e.g. Reddy et al. 2008). At 
a redshift z = 1.433, this corresponds to an observed wavelength 
of 4136 A, which falls between the transmission of the SDSS u 
and g filters. From Table [T] we find that u(il + i2) = 20.9 and 
3(11 + 12) = 20.87. Adopting a magnitude mi7oo(il+12) = 
20.9 + 2.51og(l + z) = 21.9 (on the AB scale), we deduce 
that at z — 1.433 this corresponds to an absolute magnitude 
AfiToo (11 + 12) — —23.2. Correcting for the magnification by 
a factor of 2.5, we find that the source luminosity at 1700 A is 
Li7oo — 61/*, compared to -L1700 = —20.3 obtained by inter- 
polating between the luminosity functions of star-forming galaxies 
at z = 1 and 2 (Reddy & Steldel 2009). 

Such a high luminosity may indicate that our lensing model 
underestimates the magnification of the source. However, we also 
note that in the typical z ~ 2 galaxy (equivalent data are not yet 
available for galaxies at z ~ 1.5) the UV continuum at 1700 A is 
dimmed by a factor of ~ 4 — 5 by dust absorption (Reddy & Steldel 
2004; Reddy et al. 2006; Erb et al. 2006c). As we shall see below, 
the source in CSWA 20 is unusual in showing no evidence for dust 
absorption. Thus, its high UV luminosity may well be the result of 
unusually low reddening, rather than an exceptionally large number 
of OB stars. 
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Table 3. Emission Lines identified in Images il and i2 of 
CSWA20 



Line 






Comments 


[Nil] 


6585.27 


^ 0.45 


This line is undetected 


Ho 


6564.614 


4.9 ±0.1 


Affected by sky residuals 


[Olll] 


5008.239 


7.0 ±0.1 




[Olll] 


4960.295 


2.23 ±0.07 




H/3 


4862.721 


1.65 ±0.04 




[Oil] 


3729.86 


2.09 ±0.03 


Blended with [O II] A3727 


[On] 


3727.10 


1.56 ±0.03 


Blended with [O II] A3729 


Mgll 


2803.5324 


0.27 ±0.05 




Mgll 


2796.3553 


0.48 ±0.06 




Cm] 


1908.734 


(o.isr 


Noisy 


[C III] 


1906.683 


(0.25f 


Noisy 



^ Vacuum wavelengths. 

^ Integrated line flux in units of 10^^*^ erg s^^ cm~^. 
Uncertain measurement. 



At 2: = 1.433 our data cover the rest-frame wavelengtli inter- 
val from ~ 1350A to ~ 9050 A. Table|3]lists the emission lines 
identified in the composite X-shooter spectrum of images il and i2 
of CSWA20; the most prominent of these are reproduced in Fig- 
ure |4] The presence of a number of other spectral features can be 
surmised from the data, although their S/N is too low (with the rel- 
atively short integration time devoted to this object during instru- 
ment commissioning) to warrant their measurement. Among these, 
we recognise the P Cygni profile of ClvA1549.1 due to massive 
stars, which is a common feature of the integrated spectra of star- 
forming galaxies (e.g. Schwartz et al. 2006; Quider et al. 2010). 

As can be realised from inspection of Figure|4l the line profiles 
are asymmetric with an extended blue wing, suggesting that they 
consist of more than one component. Inspection of the 2D images 
confirms the presence of two components, separated by ~ 0.4" 
along the slit and by ~ 100km s~^ in the spectral direction. It is 
possible that the source consists of two merging clumps of gas and 
stars. In future, it should be possible with better data to extract these 
two components separately and compare their properties. Given the 
limited S/N ratio of the current observations, we opted for a single 
extraction (Section |2.3l l which blends together the light from the 
two clumps, and results in the asymmetric line profiles evident in 
Figure |4] When these emission lines are analysed with Gaussian 
fitting routines, as explained below, they appear to consist of two 
components with the parameters listed in Table |4] However, it is 
important to keep in mind that these parameters refer to the blend 
of the emission lines from the two clumps, and may well turn out to 
be different from the individual values appropriate to each clump 
when the two are analysed separately. However, the total flux val- 
ues deduced for the emission lines and listed in Table |3] should not 
be affected. 

Thus, in order to deduce redshifts z^m, velocity dispersions 
a, and line fluxes F, we fitted the emission lines with two Gaus- 
sian components, using ELF (Emission Line Fitting) routines in the 
Starlink DIPSO data analysis package (Howarth et al. 2004), as well 
as custom-built software. The fitting proceeded as follows. 

The best observed emission lines among those listed in Table|3] 
are the [O III] doublet and H/3, as they are relatively strong and free 
from sky residuals. Consequently, we first fitted these three nebular 
lines, using custom-built Gaussian fitting routines which allowed 
the redshift, velocity dispersion, and flux in each component to vary 
but with the constraint that the redshift and velocity dispersion of 



Table 4. Profile Decomposition of [O hi] and H/3 Emission Lines 



Component 




(T^ (km s^i) 


1 


1.43308 ± 0.00004 


85 ±3 


2 


1.43354 ± 0.00001 


26 ±3 



^ Corrected for the instrumental resolution. 



each component should be the same for all three lines, and therefore 
using the data from all three lines to find the values of Zcm and 
a which minimise the difference between computed and observed 
profiles. Errors on the values of Zcm, cr, and F so deduced were 
estimated using a IVIonte Carlo approach, whereby the best fitting 
computed profile was perturbed with a random realisation of the 
error spectrum and refitted. The process was repeated 100 times and 
the error in each quantity (5zcm, 5F) taken to be the standard 
deviation of the values generated by the 100 IVIonte Carlo runs. 

Table |4] lists the values of 2cm, and a so derived (after sub- 
tracting in quadrature the value of ctinstr corresponding to the in- 
strumental resolution of the appropriate X-shooter spectrum — see 
Section |2l2l l. The theoretical profiles computed with the parameters 
in Table|4]and the line fluxes in Tableware shown with continuous 
lines in Figure|4] 

In the next stage of the process we applied the values of 2cm 
and o deduced from the analysis of the [O III] and H/3 lines to the 
[O II] doublet which is only partially resolved, leaving the flux in 
each member of the doublet as the only free parameter. As can be 
seen from Figure|4] the model parameters in Table |4] provide a sat- 
isfactory fit to the [O II] doublet. The fit to Ha seemingly required 
a small shift in the redshift of the narrower component, correspond- 
ing to a velocity difference Av = +25 km s~^, or ^ 2 wavelength 
bins, but we suspect that this is an artifact of sky residuals which 
affect the Ha emission line more than other spectral features. The 
integrated flux in the line is the same independently of whether this 
shift is applied or not. 

The measurements collected in Tables |3] and |4] allow us to 
determine a range of physical properties for the lensed galaxy in 
CSWR20. Before discussing the more involved derivations, we 
point out two straightforward conclusions. First, we note that the 
ratio F(Ha)/F(H/3) = 2.95 ± 0.1 is as expected from Case B re- 
combination, F(Ha)/F(H/3) = 2.86 (Brocklehurst 1971), with the 
'standard' parameters of temperature T = 10 000 K and electron 
density n(e) = 100 cm~^ (the density dependence is in any case 
minimal, and the temperature dependence is minor). Thus it ap- 
pears that the emission line gas is essentially unreddened. A lack 
of dust in this galaxy (as viewed from Earth) is further indicated 
by the blue UV continuum: the photometry in Table [T] implies a 
UV spectrum which is flat in fv, as expected for an unreddened 
population of OB stars (e.g. Bruzual & Chariot 2003)0 

Second, the [O ll] doublet ratio, which is sensitive to density 
(Osterbrock 1989), is found to be close to the low density limit: 
F(3729)/F(3727) = 1.33 implies n(e) = 110 cm"^ While at 
face value the ratio of the [C III] doublet lines would suggest much 
higher electron densities (> 1000 cm^^), these lines are consider- 
ably weaker than [O II] and their ratio is much more uncertain, so 
that we consider the value of n(e) deduced from [O ll] to be the 
more reliable. 

^ Galactic extinction is negligible in this direction, with E(B — V)mw = 
0.022 (Schlegel, Finkbeiner, & Davis 1998). 
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5.2 Star-Formation Rate 

The H/3 flux given in Table[3]implies a luminosity L(H/3) = 2.1 x 
10*^ erg in our cosmology. Adopting the Case B recombination 
ratic0, F(Ha)/F(H/3) = 2.86, the corresponding Ha luminosity is 
L(Ha) = 6. Ox 10'*^ erg s~^, which in turn implies astar formation 
rate: 

SFR = 7.9 X 10-*'L(Ha) ■ 1^.1^.2 = 21(Mq yr"!) . (2) 

The first term on the right-hand side of eq. {2) is the conversion 
factor between L(Ha) and SFR proposed by Kennicutt (1998), to 
which we apply three corrections, as follows. The first adjustment, 
by the factor of 1/1.8, takes into account the flattening of the stellar 
initial mass function (IMF) for masses below lil/© (Chabrier 2003) 
compared to the single power law of the Salpeter IMF assumed by 
Kennicutt (1998). The second correction factor is the 2.5 x magni- 
fication we estimated for the sum of images il and 12 (Section^. 
The last term corrects for light loss through the spectrograph slit, 
which we estimate to be a factor of ~ 2 by comparing the mea- 
sured UV flux in the X-shooter spectrum with the SDSS magni- 
tudes given in Table[T] A factor of ~ 2 slit loss is typical of near-IR 
observations of nebular emission lines from high-z galaxies (e.g. 
Erb et al. 2006c). 

An independent measure of the SFR is provided by the UV 
continuum from OB stars. From + 12) = g{\l + 12) = 20.9 
(Table[B, we have /^(1700) = 1.6 x 10"^* erg s"^ cm"^ Hz'^ 
from the definition of AB magnitudes in terms of The corre- 
sponding luminosity 1/^(1700) — 8.3 x 10^^ erg s"^ Hz^^ in turn 
implies: 

SFR = 1.4 X 10-''L.(1700) ■ -L ■ -L = 26(Mq yr'^) , (3) 

i.O Z.D 

using Kennicutt's (1998) scaling between Luv and SFR and ap- 
plying the same corrections as above for the Chabrier (2003) IMF 
and magnification factor. The good agreement between the SFR es- 
timates from the UV continuum and the Balmer lines is a further 
indication that the OB stars and H II regions of this galaxy suffer 
very little reddening from dust. 



5.3 Metallicity 

Since we detect emission lines of [O II], [O III], and H/9, we can use 
the i?23 = [F(3726) + F(3729) + F(4959) + F(5007)]/i^(H/3) 
index first proposed by Paget et al. (1979) as an approximate mea- 
sure of the oxygen abundance in the absence of direct tempera- 
ture diagnostics. Over the thirty years since the seminal paper by 
Paget and collaborators, many studies have shown that the values 
of (O/H) so deduced are accurate to within a factor of ~ 2 when ap- 
plied to the integrated spectra of galaxies (e.g. Pettini 2006; Kewley 
& Ellison 2008 and references therein). 

Figure[5]shows the dependence of (O/H) on R23 using the an- 
alytical expressions by McGaugh (1991) as given by Kobulnicky et 
al. (1999); these formulae express O/H in terms of R23 and the ion- 
ization index O32 = [i^(4959) + F(5007)]/[F(3726)-H 7^(3729)]. 



^ We go through this route, rather than using the Ha flux directly, because 
the Ho emission line is contaminated by strong sky line residuals and we 
consider the measurement of its flux less reliable than that of H/3. How- 
ever, in practice we would have obtained essentially the same result had 
we adopted the value of _F(Ha) from Table [3] as our starting point in the 
following calculation. 




log(Rj = log(([0 Il] + [0 III])/H,S) 

Figure 5. Oxygen abundance from the R23 index. The continuous lines 
are the calibration byMcGaugh (1991) for the measured ionization index 
O32 = [O III]/[OlI], while the shaded area shows the values allowed by the 
measured R2Z line ratio and its statistical Icr eiTor. The broken horizontal 
line is the reference solar oxygen abundance 12 4- log (O/H) = 8.69 
(Asplund et al. 2009). 

The double-valued nature of the _R23 index translates to an uncer- 
tainty in the oxygen abundance between 12 + (O/H) = 8.09 
and 8.52, or between ~ 1/4 and ~ 3/4 of the solar abundance 
12 + (0/H)q = 8.69 (Asplund et al. 2009). 

The degeneracy can be broken by considering the [N II]/Hq 
ratio. According to Pettini & Paget (2004), our upper limit on the 
N2 = \og[F{6583)/F{Ha)] index, N2 < -1.04 (Table[3), im- 
plies 12 + (O/H) < 8.26, or ^ 2/5(0/H)q, favouring the lower 
branch solution of the ii23 method. A metallicity of ~ 1/4 solar is 
not unusual for star-forming galaxies at z = 2-3 (Erb et al. 2006a; 
Maiolino et al. 2008). 

5.4 Mg II Emission 

We conclude this Section by commenting briefly on the presence 
of narrow Mgll AA2796, 2803 lines among the nebular emission 
seen in CSWA20. These lines are weak (see Figure |6ll, although 
undoubtedly real (8cr and 5a detections respectively — see Table|3]l. 
To our knowledge, narrow Mg II emission has rarely been reported 
in nearby extragalactic H II regions, but this could be due to the 
fact that this wavelength region has not been observed extensively 
in extragalactic nebulae and, if the two lines are generally weak, 
existing observations may not have the sensitivity and resolution to 
detect them. 

Broad Mgll emission is of course a common feature in the 
spectra of Active Galactic Nuclei (AGN), but we find no evidence 
in our spectrum of CSWA 20 for the presence of an AGN. Specifi- 
cally: (a) the Mg II lines are narrow (the profile decomposition into 
two Gaussians returns values of a comparable to those listed in 
Table |4); (b) we detect no high ionisation emission lines due to 
an AGN in the rest-frame spectral range 1350-9050 A of our data; 
and (c) the nebular emission line ratios fall well away from the lo- 
cus occupied by AGN in diagnostic diagrams such as that shown in 
Figure|7] 

It is possible that in star-forming galaxies at intermediate red- 
shifts weak Mg II emission is more common than anticipated. In 
their survey of 1406 galaxies at z ~ 1.4, Weiner et al. (2009) iden- 
tified a small proportion of galaxies (50/1406, or ~ 3.5%) with 
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Mg II AX2796, 2803 



X 




_ 6800 6850 

Wavelength (A) 

Figure 6. Mg II AA2796, 2803 emission lines in the CSWA20 lensed 
galaxy. The black histogram shows the data, while the blue continuous line 
is our two component Gaussian fit (see text for further details). 



'excess' Mgll emission whose nature they were unable to estab- 
lish with certainty. However, even their stacked spectrum of the 
remaining 1356 galaxies does show weak Mgll P Cygni profiles, 
with broad blueshifted absorption and weak emission redshifted by 
a few 10s of km s^^. Galaxies with 'excess' emission tend to be 
bluer than average. 

CSWA 20 appears to fit into this general pattern, with the ex- 
ception that no strong absorption component is evident. Thus, it is 
possible that the unusual lack of absorbing material in front of the 
stars and H II regions suggested by the negligible reddening noted 
above gives an unimpeded view of the intrinsic Mg II emission from 
the H II regions of this galaxy. It is interesting that, when fitted with 
two Gaussian components, as described in Section [5TT] for the other 
nebular lines, the best fitting values of redshift for the two com- 
ponents in Mgll differ by Ai; = +30 km s^^ (~ 3 wavelength 
bins) from the corresponding values listed in Table |4l We do not 
believe that this difference is due to an incorrect wavelength cali- 
bration, because we checked the accuracy of the wavelength scale 
in this region of the spectrum by reference to sky emission lines, 
and found it to be ±2 km s^^ (~ 1/5 of a wavelength bin). If the 
shift to longer wavelengths is not due to noise, it may be an indi- 
cation of radiation transfer effects in an expanding medium, analo- 
gous to those commonly seen in the Lyman a line (e.g. Verhamme, 
Schaerer, & Maselli 2006; Steidel et al. 2010). 



6 SUMMARY AND CONCLUSIONS 

We have presented X-shooter observations, which are among the 
first obtained with this new VLT instrument, confirming the grav- 
itational lens nature of CSWA 20. This system, originally identi- 
fied as a candidate from the CASSOWARY search of SDSS im- 
ages, is found to consist of a luminous red galaxy at 2:abs = 0.741 
which magnifies the light from a background star-forming galaxy at 
2em = 1.433 into four images of approximately equal brightness. 
With a velocity dispersion aions — 500 km s~^, the lensing galaxy 
is among the most massive known; the mass A/ ~ 4 x 10^^ en- 
closed within the Einstein radius (-^ 21 kpc) is responsible for the 
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Figure 7. [O III]/H/3 vs. [N II]/Ha diagnostic diagram. The location of the 
lensed galaxy in CSWA 20 is shown by the red square and left-pointing 
arrow. The small green circles are galaxies from the KISS survey (Salzer et 
al. 2005) and the small black dots are local starburst galaxies from Kewley 
et al. (2001). The dashed line shows the locus of points which Kewley et 
al. (2001) consider to be the theoretical limit for starbursts, in the sense that 
galaxies without an AGN component should fall below and to the left of 
this line. The dotted line is an empirical determination by Kauffmann et 
al. (2003) of the same limit. The presence of an AGN in CSWA 20 seems 
unlikely, given the weakness of the [N II] emission. 



wide separations (~ 6") between the four gravitationally lensed 
images of the source. 

The source blue colours satisfy the 'BM' photometric selec- 
tion criteria of Steidel et al. (2004), which isolate galaxies in the 
redshift interval 1.4 ^ z S:, 2.0. BM galaxies have been relatively 
little studied up to now. The source in CSWA 20 shares many of the 
properties of the much better characterised sample of 'BX' galaxies 
at 2: ~ 2. Its star-formation rate, SFR ~ 25Mq yr~^, metallicity 
(0/H) ~ 0.25 (0/H)g, and velocity dispersion (equivalent to two 
components separated by Av = 57km s~^, with ai = 85km s^^ 
and (T2 = 26 km s~^), are all within the range of values appropri- 
ate to star-forming galaxies at z — 2-3. The galaxy in CSWA 20 
is, however, unusual in its low degree of reddening; only ~ 5% 
of the BX galaxies studied by Erb et al. (2006b) have values of 
E{B — 1^) ~ 0, as found for this source. Evidently, we view the 
starburst from a direction along which most of the foreground in- 
terstellar material in the galaxy has been evacuated. The lack of 
dust extinction probably contributes to the unusually high UV lu- 
minosity of this galaxy, L ° ~ 6L* « , and the absence of 

b J' 1700A 1700A 

absorbing gas in front of the H II regions may be the reason why 
we can detect weak Mg II resonance lines in emission. 

The observations presented here are a good demonstration 
of the new opportunities for studies of high-2; galaxies afforded 
by the availability of X-shooter. Data secured with less than two 
hours of integration on two of the four images (with magnitudes 
g ~ r ~ 21.5) have been sufficient to characterise many of the 
salient properties of both lens and source, thanks to the high effi- 
ciency and wide wavelength coverage of this unique instrument. 
Longer observations, achieving higher S/N ratios, will in future 
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shed further hght on the nature of both the massive foreground 
LRG and the background star-forming galaxy. The gravitational 
magnification will make it possible to study the internal structure 
of the latter on unusually small scales, by comparing the individ- 
ual spectra of the four images (e.g. Stark et al. 2008). With more 
strongly-lensed systems continuing to be found, the characteristics 
of X-shooter will veritably open a new window on the high-redshift 
universe, well ahead of the era of 30-1- m telescopes. 
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